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ABSTRACT 


The Rocket Propulsion and Combustion Laboratory at the 





Naval Postgraduate School was evaluated to determine if the 


currently installed systems could support the operation of 





a subscale hypersonic test stand. Trait analysis of 
existing lab support sub-systems was performed and an 
envelope of operation, defined by flight altitude and 
flight Mach number was determined. Calculations revealed 


that if the air supply volume was doubled to 160 ft?, a 





hypersonic blow-down facility could generate conditions to 





perform testing in the dynamic pressure window of 500 to 
2000 psf, from Mach 5 to 7.5 with a maximum system mass 
flow rate variation from 3 to 45 lbm/s. Additionally, a 
dynamic design process was outlined to assist other 
designers in producing similar test stands. Finally, a 
software analysis package was developed to analyze proposed 


changes in the support system architecture. 
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pe INTRODUCTION 


A. MOTIVATION 





In August of 2002, an Australian led design team 
tested a low-cost hypersonic engine. The air-breathing 
scramjet (supersonic combustion ramjet) engine was 
accelerated to supersonic speeds by means of a solid rocket 
motor powered missile. The air-breathing engine was then 
started and positive thrust was achieved. This modest 


success has motivated designers to reinvestigate 








applications for scramjet engines. One area of interest to 





the space industry is combined cycle rocket engines, that 
promise to reduce launch costs from $10,000 per lb to $1000 
per lb. To support the verification of new designs, 
hypersonic test facilities are necessary. A forward 


thinking method of developing an optimized hypersonic test 





facility was formulated. Specifically, the Rocket 


Propulsion and Combustion Laboratory (RPCL) at the Naval 





Postgraduate School was evaluated to determine if the site 
and support sub-systems were suitable to support a 
hypersonic test stand. 

Scramjet propulsion provides an air breathing method 
of achieving speeds above Mach 5, and provides the promise 
of reducing the cost to access space. Traditional vertical 


launch methods are costly, estimated at $10,000 per pound 





of payload because the fuel and oxidizer must be carried 


onboard the vehicle, and accelerated, with the payload to 





escape velocities. Additionally, bipropellant rockets 








provide limited specific impulse (Isp). A combined cycle 
launch vehicle may contain a hypersonic engine as an 


intermediate stage on a conventional chemical rocket. Such 


a combined cycle configuration would use the chemical fuel 


and oxidizer stored on board the vehicle for initial launch 








and acceleration to speeds where sustained supersonic 
combustion can be achieved at, approximately Mach 4 to 5. 


At this time, the air breathing engine would be used to 





augment the thrust of the main rocket engines. This 








concept provides a higher Isp and a lower cost to access 
space due to the fact that the combined cycle rocket 


carries less fuel and considerably less oxidizer and thus 





its mass is significantly less than a conventional launch 





vehicle. This drastic reduction of the launch vehicle’s 





mass yields an order of magnitude increase in payload mass, 


which reduces the cost to launch payloads into space down 





to roughly $1,000 per pound (Ref. 1). 





Hence, the desire to test small scale hypersonic 


engines for sustained periods of time provides the 





motivation for universities to establish low-cost test 








facilities to study the processes involved in optimizing 


the performance of thes ngines. This thesis provides a 





road map to the process of designing and optimizing a blow- 
down facility for the testing of subscale combustors. Air 
heating methods are introduced but, vitiated heating is 
analyzed for this design, due to its simplicity and low 


cost. 


B. BACKGROUND 


There is a wide variety of methods available to 





achieve the supersonic flows required to supply scramjets, 





and as Lu and Marren (Ref. 2) state, that the goal of any 


ground based test facility is “Duplication” in which all 





aspects of the flight environment are matched. Lu and 





Marren note that this is rarely achievable and further 
2 


state that useful test results can be derived from a 





facility that “Replicates” the flight environment. In this 


realm, flight temperature, pressure, velocity and chemical 





composition conditions are met. Finally, “Simulation” is 
presented as the last area of operation where some 


important parameters of the flight environment are equated. 


1. Air Heating Methods 





The cost of the test facility is directly dependent on 








the accuracy to which the flight environment is duplicated. 





For example, chemically accurate air compositions are most 
readily obtained by heating the air with a pebble-bed 


heater, but these units are costly and the university test 








facility designer is forced to utilize other means. iG. <2 








important to note that achieving the high enthalpy required 





to match the flight environment is perhaps the most 





challenging aspect of designing a hypersonic test stand. 
The design of the heater affects the maximum enthalpy of 


the airflow, as well as the chemical composition of the air 





and the maximum mass flow rate sustainable by the system. 


Since so many test parameters rely on the design of 
this component, it is wise to investigate some of the 


available options further. 


a. Pebble Bed Heating 





As stated earlier, one option is th pebble-bed 
heater, which consists of a several layers of ceramic 
spheres that are heated using a combustion cycle or 


electrical heating units. When a target temperature is 








reached, the combustor is turned off and high pressure air 


flow is directed through the matrix of pebbles, where it is 


3 








heated. In this manner, no chemical impurities are 
introduced to the air and the chemical composition is 


accurate to what would be seen during in-flight conditions. 








However, as stated before, these units are prohibitively 


costly for most universities. 


b. Vitiated Heating 
Combustion heating of the supply air is another 


option that is available to produce the high enthalpies 





needed to support testing. In this method, a high pressure 








burner section located before the facility nozzle is used 


to start, support and maintain combustion during the 





duration of the test run. The chemical composition of the 


air supplied to the facility nozzle will depend on the type 





of fuel used to heat the air. When hydrogen is used as the 
fuel, with oxygen and N2 as the oxidizer, water and very 
minute amounts of NO,, OH, Os, No, He, O, and H will be the 


by-products of the combustion. The concentration of each 





of thes species depends on the conditions in which the 





combustion occurs. The result is essentially air with by- 








products and surplus water which absorbs a significant 
amount of heat. Vitiated heating as this method is called 


also requires running the fuel and oxidizer flow-rates at 





stoichiometric correct ratios to ensure that the O») content 








in the final mixture is correct. Depending on the size of 





the facility and therefor th flow-rates involved, this 


method provides acceptable facility run times, and the most 





economical solution to the university facility designer. 


c. Arc and Shock Heating 

Two other air-heating methods are not suitable 
for this analysis because they support only very short test 
durations which is unsuitable for testing the operation of 


the combustor sections of hypersonic vehicles. 


Arc heating of the airflow is one of these 


methods and while it is capable of achieving very high 





enthalpies, Lu and Marren state that chemical impurities 








will most likely be present in the air and only small flow 





rates are sustainable. Finally, shock facilities are 





capable of producing “clean” air but the test time length 


is not sufficient for testing combustor sections alone. 





In this thesis, the use of vitiated heating was 





chosen as the method used for generating the high 





enthalpies needed to support this test facility. Vitiated 


heating reguires a lower initial investment than pebble-bed 











heating and if designed and operated correctly, it can 


produce high guality air over a range of enthalpies and at 





a flow rates that are useful for testing scramjets. 


2. The Evaluation Process 
The first step in designing the test facility is 


determining the parameters that need to be matched and the 





ranges of interest for those variables. In the case of 
hypersonic propulsion, the parameters of interest are 
outlined in reference 2. For a blow-down facility, the 


component of interest is the supersonic combustor of a 








scramjet engine, and the combustor inlet conditions are the 


parameters that need to be matched. However, depending on 





the specific scramjet design and vehicle flight profile, 





the combustor inlet conditions vary greatly between 
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designs. In an effort to ensure maximum utility of the 
final design, the collection of performance data, by Heiser 
and Pratt (Ref. 3) was consulted. Specifically, the data 


for combustor static pressure, static temperature and Mach 





number were used. 


a. Combustor Pressure Limits 





Heiser and Pratt identify the following limits. 





Combustor static pressures ar typically in the range of 











0.5 to 10 atmospheres. The upper limit is ™.. a value that 
will lead tO “reasonable” or “acceptable” weight, 
complexity and cost.” While the lower limit is to minimize 
the “.. length of the combustor required to complete the 


reaction and consume the available fuel..”. 


b. Combustor Temperature Limits 
The combustor temperature limits are based on 


maximizing the thermodynamic cycle efficiency. For flight 





at the same altitude and Mach number, a lower combustor 


inlet temperature means additional heat must be added in 








the combustor to support the same thrust. This means 
greater fuel consumption by the engine and lower 
efficiency. Therefore higher combustor inlet temperatures 





are desirable, but as Heiser and Pratt note, if the 


combustor inlet static temperature is too high, 





dissociation will occur, which results in a loss of energy 


due to the dissociation process and subsequently lower 











efficiencies. Because of this limitation, Heiser and Pratt 
state that the maximum combustor entrance Static 
temperature is ™.. almost always found to be in the 
relatively narrow range of 2600-3000 °R..” To bound the 
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value for temperatur requirement, their ™“..representative 


estimate of 2800 °R..” was used. 


c. Combustor Mach Number Limits 


Finally, appropriate values for the combustor 





inlet Mach number wer needed and determined again from 





reference (3). A close look at the relationship between 





Mach number and static temperatures shows that for air at 





temperatures ranging from ambient to 3000°R the average 
value for ?, the ratio of specific heats, is 1.36. The 


maximum ratio of combustor to free stream static 





temperatures is normally found to be about 7. Using these 
values, the minimum flight Mach number that ensures 
supersonic combustor inlet flow was calculated to _ be 


approximately Mach 7. 


C. OBJECTIVES 

This collection of data, summarized in Table (1), 
forms the boundaries within which a test facility should 
operate in order to supply flow useful for testing 


hypersonic combustors. 





Table 1. Combustor Inlet Boundary Conditions 





Based on this information, the volume of the air tanks 
and the size of the supply lines and limitations of 


installed regulators at the RPCL were analyzed to determine 








if they could support the operation of a test facility with 





these outlet conditions for a period of time that would 
allow the collection of meaningful data to evaluate sub- 
scale tests of hypersonic combustor designs. Additionally, 
this information formed the basis for a program that was 
developed to allow for analyzing the effects of changing 


many different parameters. 


II. DESIGN AND ANALYSIS 


A. PHYSICAL LIMITATIONS 


The first calculation in the process was to define the 








flight envelope; this condition limits the stress on the 





vehicle frame at higher dynamic pressures and ensures 





sufficient oxygen concentration at lower dynamic pressures. 


The values were calculated using equations (1) through (4) 





and the United States Standard Atmosphere Tables 1976 (Ref. 
4). 


us 
a 


a= ¥g.RT (2) 


P=pRT (3) 
pv’ 
oo 4 
dynamic 2g. ( ) 
Ph = T PagicM : (5) 


where: Pgstatic 1S a function of altitude 


The formula for dynamic pressure is equation (4) but 





the density is altitude dependent and the velocity is 





temperature dependent. To account for these variations, 








equation (3) was solved for density and equation (1) and 
(2) were solved for velocity. Resulting in equation (5), 


the dynamic pressure was then calculated as a function of 





the independent variables which were the Mach number 


defined in eguation (1) and Petatic. The next aspect of the 








solution was to determine the combustor stagnation pressure 
9 


and temperature conditions that fall inside this flight 








profile. To do this the stagnation conditions were 
calculated using the fr stream conditions using equations 
(6) and (7). Then those results were translated using 





equation (8) and the fact that the adiabatic compression 





process does not change stagnation temperature conditions. 











-] -1 
P= Pay {+r} (6) 
-1 
T= T atic 10 Lye (7) 
The fr stream stagnation temperature was calculated 





from the static temperature which was a function of 





altitude, obtained from the atmospheric model, and the free 





stream Mach number which was known. Since the compression 


process using oblique shocks was analyzed as an adiabatic 





process, the free stream stagnation temperature and the 





combustor inlet stagnation temperature had to be equal. 


Also one of the premises of the problem was that the 








maximum combustor inlet static temperature was limited to 


2800 °R, (Ref. 3) as stated previously, this conservative 








value was then substituted into the stagnation temperature 








relationship and used to find the combustor inlet Mach 


number. 


The fact that the oblique shock compression process 


is adiabatic means that the fr stream and combustor inlet 











stagnation pressures are not equal. Therefore, quation 


(8) is used to relate the static conditions in the free 
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stream and the combustor inlet through the use of the 











compression efficiency, Lees and the maximum static 
temperature ratio, ?. To perform this calculation, the 
free stream static pressure was determined next. It is a 


function of altitude and obtained via the atmospheric 
model. The free stream Mach number was the independent 
variable and therefore a known value. Next, the combustor 


static pressure was evaluated using equation (8). 


W (y-1) 
= = 8 
combustor freestream wv (1 i: Nn, ) + Nn, ( ) 





Then combustor inlet stagnation pressure was 


calculated using the combustor inlet Mach number and 





equation (6) the stagnation pressure isentropic 





relationship. 


B. EVALUATION OF RPCL’S SCRAMJET FACILITY 


With the physical limitations of the processes 





defined, the next step was to evaluate the proposed site. 


This process of designing the components to supply the high 





temperature, high pressure and subsequently the high Mach 
number air to a test combustor can be accomplished in two 


different ways. 





In the first method, the air supply requirements of a 
prototypical scramjet are supplied to a test facility. 
Normally this is in the form of stagnation pressure, 
stagnation temperature, Mach number and mass’ flow-rate 


requirements. Using this information an engineer can 
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design the heating section and convergent-divergent nozzle 


to develop these flow conditions. The result would be a 








test stand that is specifically designed to support such a 








prototype combustor. It is however unlikely that this test 


stand would fully utilize the capabilities of the facility 





support systems (i.e. specifically, the air supply system). 


Additionally, et: would probably reguire significant 





redesign to accommodate future prototype combustors that 





were submitted to the test facility for evaluation, even if 
only one of the test parameters that formulated the initial 


design were changed. 


A more dynamic approach allows the designer to define 


the operational envelope of a proposed test stand based on 





current laboratory support system capabilities. In this 


way, a proposed solution can be designed to utilize the 





full capacity of any existing facility support systems. If 


the design process is translated into software, then the 








designer gains the ability to quickly investigate the 
effects of changing certain parameters on the performance 


of the proposed test stand. Allowing this evaluation of 








performance without having to purchase additional hardware 














is a very attractive feature to designers working with 





small budgets. In the end, a suitable and versatile 








solution to testing sub-scale hypersonic combustors can be 


achieved. Using this second method, the RPCL, at the Naval 








Postgraduate School was evaluated to determine its ability 


to support the operation of a hypersonic test stand. 


1. Existing Subsystems 
Several of the required support sub-systems were 


already in place and operational. These included, the air 
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supply system, and hydrogen and oxygen supply systems. The 
air sub-system is supplied by two air compressors. Both 
are manufactured by Bauer, the compressor with the highest 


capacity is rated at 50 cfm and the other delivers 





approximately 30 cfm. This high pressure air is directed 
to four 20 ft* high pressure storage tanks which can store 
air at a maximum pressure of 3000 psi. Air from these 


storage tanks is piped to several test bays using 1%” 








stainless steel tubing. A pressure regulating valve 





upstream of the test bay connections regulates actual air 


supply pressure up to 1000 psi. Compressed hydrogen and 








oxygen cylinders are also installed at the facility and 








directed to several test bays using %” stainless steel 


tubing. 


a. Air Subsystem 


The choked air flow rate for the current air 





supply piping size was calculated to generate an upper 


limit for the maximum mass flow rate that can be supplied 





to a test article. 








It is important to note that the air tank static 


temperature changes over the course of a test run. This is 








due to the fact that the air is being supplied through an 


isentropic expansion process from the air supply tanks. 





This aspect of the process was included to account for the 
large variation in supply temperature that occurs during 


tank blow-down. 


Starting with the isentropic expansion equation 


to determine mass flow rate/area: 
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The variables above wer specified and supplied 





to the equation and a choked area for the desired mass flow 
rate was found. The process was repeated for all of the 


gas supply systems using the appropriate system parameters. 





Equation (9) was evaluated for the values of 
combustor stagnation pressure that fell within the flight 


profile. Stagnation pressures ranged from 337.5 to 2563 





psi and the stagnation temperature was set to 70 °F. This 





ensured that a test stand designed to accommodate this 
maximum mass flow rate would not encounter choked flow in 


the air supply system. The air piping diameter is 1 %” 





downstream of the regulator. A high stagnation pressure 





estimate of 2563 psi with a stagnation temperature of 70°F 
(530°R) was used to calculate the maximum air system mass 


flow rate. As stated before, the temperature estimate is 





conservative because no cooling was assumed in the choked 





flow calculation. Using these parameters, ASE was 





determined that the maximum airflow rate was about 50 


lbm/sec. 


The choked system flow rates were then calculated for 





each of the combustor stagnation pressure and temperature 





values that were found to fall within the flight window. 


The results for the total air system flow rates varied as 





shown in figure (1). 
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Mass Flow Rate vs Stagnation Pressure 
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Figure 1. Mass Flow Rate as a Function of Temperature for 2560 
psi 


b. Vitiated Heating Subsystem 


The next calculation that was performed was based 





on determining the required mass flow rate to generate the 
combustor inlet stagnation temperature at the combustor 
inlet stagnation pressure. This was done using a look-up 
table based on combustion data calculated with Thermo- 


Chemical Equilibrium Program(TEP). The combustion data was 





generated using hydrogen as the fuel and varying amounts of 
air and oxygen as the oxidizer, and was performed at 1/7, 
34, 51, 68, 85, 102, 120 atm. The air mass concentration 
was held constant and the fuel/air ratio was varied to 
produce a mixture with 21% oxygen after combustion. The 
final oxygen content was adjusted so that all variations 


were contained within +0.33%. 
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Variation in O, Content in Vitiated Air 
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Figure 2. Variation of Oxygen Content in Vitiated Air 
[Pe=1000 psi] 
The air mass concentration and fuel/air ratio were 
recorded. The mixture’s molecular weight, density, and 


ratio of specific heats were also tabulated and entered 
into the look-up table. Then using the relationship 
between mass flow rate of air and the air mass 
concentration, the mass flow rate of the oxidizer and 
oxygen were found. Similarly, using the mass flow rate of 
oxidizer and fuel/air ratio, the mass flow rate of hydrogen 
was found. These relationships are described in equations 


(10) through (14). 


m MN) xidizer 


total ~~ 


+ M fet (10) 


m =m,,+m (11) 


oxidizer ~ "" “air oxygen 


16 


Mai = Weight % ji, (ii ceeds ) (12) 


mM . 
n .. =— rr 13 
MN oxidizer Weig h t% - ( ) 


M2 ve a Miyarogen =F/ A(Tit,. con) (14) 


These values were calculated so that the hydrogen and 


oxygen flow rate limits were not exceeded. 


2. Testing Time 


With the air system flow rate maximized, it was 








necessary to determine the run time that the air flow rate 
could be supplied for. Therefore using an air system mass 


flow rate, 3000 psi initial supply air pressure and a 





starting supply temperature of 70 °F as the initial 
conditions in the air supply tank, the corresponding system 


run time was calculated. This was done using the perfect 





gas relationship, and the definition of mass flow rate. 





The mass flow rate equation was simplified for steady flow 
conditions, which is an assumption of the calculation. 


Additionally, the continuity equation and the isentropic 





expansion process were used. 


dm 
alla 15 
aaa (15) 
Am 
1 = — 16 
oe (16) 
ih = Min itial — TN final (17) 


t 
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m=pAV (18) 





Vinal = P pital (19) 
T pitiat Pinal 
= PinitiaaY Z P ina (20) 
m 


Finally the equations were solved for system run time. 
This result was then evaluated to determine if it was 
sufficient for conducting meaningful testing based on the 
criteria outlined by Pope and Goin (Ref. 5). Specifically, 


they assert that many factors should be considered when 





determining acceptable system run times. They conclude 


that the facility designer needs to address pressure 








stabilization inside the test stand, instrument response 





times for the test equipment and the type of test to be 
run. While the instrument response characteristics have 
undoubtedly improved since 1965, the rate at which pressure 
stabilizes and system test profiles are still salient 
points of consideration. Thus this thesis used their 


initial estimate for acceptable system runtime of 30 s. 





If the system run time result is less than 30 s, then 








the run time was set to 30 s and the air system mass flow 








rate is calculated, and used for further calculations. 


3. Facility Nozzle Throat Area 
The sum of the oxygen, air and hydrogen mass’ flow 


rates wer then used to calculate the facility nozzle 





choked area, using the choked flow equation (9) presented 


earlier. Since the process in the facility nozzle is 
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isentropic, the vitiator stagnation pressure and stagnation 





temperature values were used in this calculation as well, 


refer to equations (6) and (7) (From Ref. 6). 


The facility nozzle supersonic area ratio was 


calculated as follows: 

















a ‘ 10 yu Pr) 
AM (y +1) (21) 
2. 

Equation (21) was evaluated to determin th area 
ratio required to generat th highest combustor Mach 
number. Specifically, the solution was found using 
Newton’s Method (From Ref. 7). Using a starting guess, of 


Mach 5, equation (21) and (22) were evaluated. 














= 2(y-1) 
2. 
d 
M (y +1) 
A 
(4) 
= (22) 
dM dM 
Then the Mach number, equation (21) (the function), 
and equation (22) (the derivative), were substituted into 
equation (23). 
M 
M,,=M,- at (23) 
df (M) 
d(M) 


When the Mj11;=0, the solution was the M; used in 
evaluating equation (21) and (22), otherwise M;=Mii1, the new 
Mach number estimate and the process repeated until a 
solution was found. Convergence occurred in 2 iterations 
for almost every case due to the well chosen starting point 


and the properties of the function itself. 


The values for facility nozzle throat area at the 
various stagnation temperatures and the 2560 psi stagnation 
pressure (the most limiting condition) were plotted as 
shown in figure (3). The facility nozzle design was then 


based on these values. 
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Figure 3. Facility Nozzle Throat Area for 2560 psi 


C. TEST STAND DESCRIPTION 
Figure (4) is a schematic of the major components 


required in the construction of a blow-down test facility. 
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The supply tanks, vitiator, and facility nozzle were 
analyzed in the previous section. Although various heat 


transfer processes are listed in figure (4), not all of 














these were modeled. Further refinement of the presented 
analysis should consider the modeling of these additional 


processes. 


Nonetheless, all variable targeted conditions depicted 








at the bottom of figure (4) were computed in the presented 


























analysis. The typical values are listed in table (2). 
Table 2. Matched Conditions 

Property Value/Range Units 

Stagnation 337.97 2.9.60 psi 

Pressure 

Stagnation 2960-4400 °R 

Temperature 

Mass Flow Rate 3-45 lbm/s 

Run time 1-30 Ss 
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Figure 4. Process Overview (not all processes were modeled) 


The thermodynamic processes that occur are further 
defined in the enthalpy versus entropy plot shown in figure 
(5). Although the open cycle is not drawn to scale it does 
show the important processes and the legend shows the 


correlation between flight and test stand states. 
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Figure 5. Schematic Representation of the Thermodynamic 
Processes of the Test Stand (not to scale) 


D. SCALING 


Another aspect of the design is scaling. The test 
article will be some fraction of the size of the flight 


combustor and therefore the flow must be appropriately 





adjusted to account for the smaller physical size. The 


test article needs to be both geometrically and dynamically 





Similar to the flight article. Once the geometric 





Similarity requirement is met, dimensional analysis is 
performed to identify parameters that will ensure dynamic 


similarity. In Ref. (8) Penner identifies the three 











factors that will ensure dynamic similarity. They are the 


Reynolds number, Damkdhler’s first number and the Mach 





number. 
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Re = —— (24) 


De ae. 
D, = : (25 ) 
2 
M= Py (26) 
\ YP 


p L =constant (27) 
Prost Loss = P ich pignt (28) 


where: if Leest = % Leight 


Een. Pecan = 2) Peg 





The dynamic similarity is then established by meeting 





the requirements of equation (27), which can be expanded to 





the version shown in equation (28). Specifically, in the 











cas wher th test article is half the length of the 











flight combustor, the test flow density must be twice the 





flight density (From Ref. 9). Applying this relationship 














to equations (24) through (26), it is easily shown that the 








Reynolds number of the flow in the test article is matched 
when binary scaling is applied, also since density 


increases with pressure at a near linear rate, Mach number 





Similarity is achieved. However the effect of binary 





scaling on Damkohler’s first number is not so readily seen. 





Damkoéhler’s first number, shown in equation (25), is a 





relationship between the combustor length, flow velocity 








and chemical conversion time. The numerator represents the 





residence time, and it is simply the time that air is 


available for interactions in the combustor. The chemical 
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conversion time however is process dependent, and is a 


function of the combustion, dissociation, and recombination 








reactions taking place. It is important to note that the 





combustion and dissociation chemical conversion times are 





the important factors in the chemical conversion time, 


because for specific mixture states they are almost always 








less than the residence time of the flow in the combustor. 





Of less importance is the time required for recombination, 








Since it is significantly greater than the residence tim 











given a flight vehicle length. It then follows, that 


Damkoéhler’s first number describes the physical state of 





the combustion reaction. 


The chemical conversion time is inversely proportional 





to a combination of the previously stated reaction rates, 





and the chemical reaction rate is directly proportional to 
the pressure at which the reaction is taking place. Since, 
binary scaling is achieved by increasing the static 


pressure and therefore the density, the chemical reaction 








rates also increase and the chemical conversion time 





decreases. Referring again to equation (25), if the length 


of the combustor is halved, then the chemical conversion 





time must be halved as well. This means that the reaction 








rates must be doubled and therefore the pressure must also 


be doubled, which was the case to increase the density to 











meet the requirements of matching the Reynolds number and 


the Mach number. 





In summary, a test article geometrically similar to 





the flight model will be dynamically similar if the flow 








Reynolds number, Damkodhler’s first number, and the flight 


Mach number are matched. This dynamic similarity is 





achieved by using the binary scaling relationship shown in 
25 


equation (277) < 


The 


binary 


scaling factor will vary 





depending on the scale of the 


stand flow 


can 


be 


adjusted 


pressure at the facility nozzle 


density. 
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test article, and the test 


by increasing the static 





outlet, which increases the 


IIIT. SOFTWARE DESIGN 


A. SOFTWARE DESIGN OVERVIEW 
The purpose of the software is to allow the facility 


designer to have a resource that will perform much of the 





comparisons that are required to determine the 





interdependence of parameters. The software was developed 





using MATLAB student edition version 6.0 and SIMULINK 




















version 4. SIMULINK was used to model most of the 





processes which allows the user to more readily see the 


flow of the solution. The MATLAB code simply links the 











SIMULINK models. Figure (6) shows the program flow path. 
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‘Combustionflow’ rate. 
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Figure 6. Program Flow Path 
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The flight conditions in the form of dynamic pressure 
and combustor stagnation conditions are calculated to 
establish the boundary conditions that the vitiator and 


facility nozzle were designed to using the 





‘atmosphericsolver’, ‘stagnationpressure’, ‘stagnationtemp’ 


and ‘Combustionflow’ models. Then the limitations of the 





installed support subsystems at RPCL are weighed into the 


design considerations in the form of sustainable mass flow 








rates and the test time using the ‘systemruntime’ and 
‘makesystemruntime’ models. Subsequently, the combustor 
stagnation conditions, combustor Mach number and mass flow 


rates are used to determine the facility nozzle choked area 








and supersonic area ratio of the facility nozzle in the 


main program and ‘findmyexitmachnumber’ model. 


B. SOFTWARE LIMITATIONS 





As with any model, limitations will be present. In 





this software, there are two categories where inaccuracies 








will arise. The two areas are physical modeling errors, 


and programming limitations. 





1. Modeling Limits 
When modeling a system, the designer has to balance 
the level of fidelity (necessary to provide meaningful 


results) and the level of effort or time to achieve them. 











Despite th ffort to ensure a high fidelity model, the use 


of one-dimensional isentropic flow is not entirely 











accurate. Although, Computational Fluid Dynamics would 





provide a more realistic description there are a few 





reasons why this path was not considered. First, the 


amount of computational time required would increase by at 
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least 3 orders of magnitude. This alone would have limited 


the number of conditions presently considered. Secondly, 





CFD requires a detailed geometric description of the 





proposed facility. This level of detail requires a priori 





knowledge that can only be provided by the models under 


consideration in this thesis. Therefore, once the final 





design is reached with the present models, one can use CFD 





to refin th design. Hence, macroscopic performance 


results were determined to be the best way to model this 





system. 


Another limitation present in this software was that 


the mass capacity of the hydrogen and oxygen supply systems 








was not considered. It is therefore possible that the 
software would report a feasible solution, when actually 
the condition is beyond the available total mass of oxygen 
or hydrogen. Accounting for this aspect of the design was 


not included in the software because of the objective to 





keep the software simple in its design. Additionally, for 








small facilities, the addition of extra O2 or H,z cylinders 
is a trivial matter which is overcome simply by connecting 


additional cylinders. To estimate the mass of oxygen 





required, the user only needs to multiply the applicable 


flow rate and the final test time. 


Additionally, the data used to describe the real gas 


effects on different air and gas mixture properties was 





organized into look up tables. The tables use linear 





interpolation for in-between values and linear 











extrapolation for values outside of the table limits. 


While the data is almost entirely linear for temperature 








and pressure values inside the table, dissociation effects 











are seen near the upper temperature and lower pressure data 
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point of the table. These values are 4400 °R and 337.5 psi 
respectively. Thus extrapolation of data values beyond 


these conditions should be viewed with caution. 


2. Programming Limitations 
Programming inconsistencies can introduce some errors 
as well and one method utilized frequently in this software 


is the use of the “hit/crossing’ function in some of the 














SIMULINK models (From Ref. 10). This function is used to 


identify the point where the ramped input, such as the 





altitude model depicted in figure (8), provides a solution 








to the problem; however, if the step size of the input 


variable is not sufficiently small, then the solution may 





actually not be the correct value. 


Regardless of the errors involved due to these 
limitations, the results of this type of analysis are 
sufficient to make preliminary design decisions, and 
highlight areas where a more detailed analysis may be 


required. 


Cc. MODELS 
1. Atmospheric Model 


When the program executes, the first figure displayed 





on the screen is calculated using the ‘atmosphericsolver’ 











model. This model generates an Altitude versus Mach number 








graph based on the limitations outlined in Ref. (2). Look- 





up tables are used to provide the state conditions of the 
atmosphere and the data used in the temperature and density 
look up tables was obtained using the United States 


Standard Atmosphere 1976 tables (From Ref. 4). The values 





in the look-up tables were compared with the 1962 United 
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States Standard Atmosphere tables included in Ref. (11) to 
ensure accuracy. The look-up tables contain data for a 
height of 60 km, or nearly 200,000 feet. But as seen in 
the dynamic pressure formulation, the flight window defined 


by Ref. (3) peaks out at about 130,000 ft. Therefore, to 





reduce the program run time, 130,000 ft was chosen as the 
maximum altitude used for calculations. The model 
subtracts a reference dynamic pressure (From Ref. 3), from 
the calculated dynamic pressure and keeps track of that 
difference. When that difference equals zero (actually 
crosses zero) the model stops and the altitude is reported, 
the program then continues execution for the next Mach 


number. Th referenc dynamic pressure values (500 and 





2000 psf) are sequentially passed to the 














‘atmosphericsolver’ SIMULINK model depicted in figure (8). 





This model is located inside of a loop that increments the 


free stream Mach number to the user defined maximum value. 











For each Mach number value, the model ramps the altitude 
input from 0 to the maximum altitude (also a user defined 
value.) Subsequently, the line of constant dynamic 
pressure is plotted on an Altitude vs Mach number plot as 
shown in figure (7). The loop continues for each reference 


value of constant dynamic pressure that was defined. 





Figure (7) shows the plot of lines of constant dynamic 
pressure. The upper line represents 500 psf and the lower 


2000 psf. 
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2. Combustor Stagnation Pressure 





The flight combustor stagnation pressure is determined 
as depicted in figure (9). Two nested loops allow the 


reference values of stagnation pressure to be varied 





between the user defined values, the upper limit being set 


by the maximum pressure in the supply air tanks at the test 





facility. The inner loop is the flight Mach number and is 
varied as described in section (a) for the 


‘atmosphericsolver’ model. The ‘stagnationpressure’ model 





shown in figure (9) uses these values to determin th 





calculated combustor stagnation pressure. 





As before, the flight altitude variable is a ramped 








input and the model uses the same atmospheric model. The 


static pressure look-up table is used to find the static 





pressure at the given flight altitude. Using equations (6) 
and (8), the combustor entrance Mach number and _ the 
combustor entry static pressur ar determined. The 








combustor stagnation pressure is then determined and the 





reference stagnation pressure (convenient values specified 
by the operation of the outer loop, and chosen within those 


bounds) is subtracted from it. When the difference equals 





zero, the model stops and control is shifted back to the m- 





file. Additionally, the altitude, atmospheric density, 
atmospheric temperature and combustor Mach number are also 


passed to the MATLAB m-file. Next, the ‘atmosphericsolver’ 





model is used again to determine if the altitude found by 





the ‘stagnationpressure’ model falls within the flight 





profile window. If it does, a point is plotted and the 





stagnation pressure is recorded. 


35 


voHeNWIS 
doys 


zadoos, 


Bussola 
HH 


cy] 
a 


gadoos 


{uonipuog 10yEHIAy 


[isd]ainssaud 
uoneubeys 
aouslajey 


Usd] sinsseid 
uoneubeys 
1oysnquioa 


(ny 


[woooos-a) #1984 dno) [4] 
ainjesadwa} ouadysoupy 
Ing fo sinjesadway [y] 0} peauog 


<-——_ 


[w go009-0) #19&} dn-yoo) 
[eyvu/B>] Qisuap ing 
ing 40 fqisuag [eaywai] op panueg 








—<: 





Jequinu yoryy 





aoedsmonn OL 7814] JO}sNquiog 
(ny 
sequinu yoeyy yY6I14 
One! sinjeradwiay 
[ols By dH ul eyeg Wold] aneys ajay 
£20 = uolssaidwios | 
£90 =suoissaidwos Z 
89'0 =suolssaidwos ¢ 
ne ewwe Z6'O =suoissaidwos p 
(Spoys enbijgo yo saquinu) 
[sL°¢ uonenba] Qua uoissaidwog ap1yay, 
[isd] ainsseud 
onEIS 
Joysnquiog 
(ny 


[w aooog-al 2192} dn4oo) 
[isd] 0} Usd] 07 [e d]ainsseid suadysowpy 
paaueg = anueg 


a 


Aeidsig apne 


9¢"L SI syeay oyloads Jo oney 
ainsseid uoeuBbe}s JUe}SUOD UO Saul] Sayejnayeo sIyL 
YOLVINITVI AYNSSsdd NOILYVNOVLS YOLSNANOD 





[w] op peau0g 


< 


Gy) spray 


SPN 














SIMULINK model 


‘stagnationpressure’ 


os 


Figure 


36 


The stagnation pressure is recorded so that it can be 


evaluated later in the program to see if it meets the 





minimum run time requirement. Also, the combustor Mach 
number is recorded if its value is larger than any previous 
values of combustor Mach number. This is done to define an 
iteration limit, used later in the program for calculating 
the supersonic area ratio of the facility nozzle. This 


process is repeated until all values between the stagnation 





pressure limits have been evaluated for each Mach number 
specified by the inner loop. These results are shown in 
figure (11). For clarity, an unreduced set of data spaced 


at 250 psi intervals is plotted in figure (10). 
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Flight Envelope 





3. Stagnation Temperature 


The next step is to determine the combustion 





stagnation temperatures that fall within the flight profile 
dynamic pressure limits. This is performed by the 
‘stagnationtemp’ model shown in figure (12) and the results 


are shown below in figure (14). 
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‘stagnationtemp’ 


V2.3 


Figure 


The ‘stagnationtemp’ model differs from the two models 





already considered because in this case the flight Mach 
number is the ramped input variable. This is due to the 
fact that the atmospheric temperature profile as a function 


of height has multiple altitude solutions for a given 





temperature value, as shown in figure (13). 
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Figure 13. Lines of Constant Stagnation Temperature 


This model has inner and outer iteration loops. The 
outer loop varies the stagnation temperature between the 
upper and lower limits defined at the beginning of the 


program. Also, the flight altitude is incremented from 
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zero to the altitude limit. The model performs a 








comparison between the calculated and reference values and 








terminates the iteration when a zero crossing is detected. 





The flight Mach number, combustion Mach number, atmospheric 
density, and atmospheric temperature are then output to the 
main MATLAB program. The altitude from the 


‘stagnationtemp’ model was checked against the 





‘atmosphericsolver’ solution of altitude for the dynamic 





pressures used to define the flight profile. Tf the 








altitude from the ‘stagnationtemp’ model is in the 





allowable range then a point is plotted on the figure, see 





figure (14). The combustor stagnation temperature is 


recorded and the combustor Mach number is evaluated to 





determine whether this value becomes th new limit for the 








supersonic area ratio calculation that is performed later 





in the program. The process was then repeated for all 





stagnation temperatures and for the respective altitudes up 





to the defined limit. 
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ure 14. Lines of Constant Stagnation Temperature in the 


Flight Envelope 





Combining the results for stagnation temperature and 


pressures at the combustor inlet creates a matrix of 


the region 





temperature and pressure conditions that defin 


which is shown in the figure (15) below. 
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4. Choked Flow Calculator 





Th pressur values are used to determine choked 


system flow rates for each of the supply gases whose 





properties are defined by the user. For each gas system, 





the ‘Chokedflowrate’ SIMULINK model, shown in figure (16), 











determines the mass flow rates. The gas mass flow rate is 





supplied as the ramped input variable. Using equation (9), 





the choked flow equation, and th user defined supply 
system characteristics, the model iterates through mass 
flow rates until the zero crossing is detected and then 
stops. The mass flow rate is recorded and passed to the 
main MATLAB program. The process repeats for each gas 
system and the results form the flow rates upper boundary 


conditions. 
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5. Vitiator Flow Rates 


Next, the ‘Combustionflow’ model represented in figure 











(17) is used to determine the mass flow rates required to 











generate the stagnation temperatures found in section (c). 





In this model, the pressur and temperature and the loop 





variables are applied to a series of look-up tables that 








determin th air mass concentration and Fuel/Air ratio 





needed to generat th specified temperature at the 





specified pressure. These values along with the air mass 








flow rate which is a ramped input are then used to 





calculate the required oxygen and hydrogen mass flow rates, 


used to generat th input temperature. The model will 








continue to increase the air mass flow rate until the air, 
hydrogen or oxygen mass flow rate reaches 80% of their 
respective choked flow values, which were found by the 


‘Chokedflowrate’ model. 


When the ‘Combustionflow’ model finishes all 
calculations, the final air, hydrogen and oxygen flow rates 
as well as the mixture gas constant (Rmx), mixture density 


(?mix) and the ratio of specific heats for the mixture (?mix) 








are reported to the main program. Tf the ©O2 or Hz mass flow 





rates were the limiting factor then the corresponding air 
flow rate is assigned as the maximum air flow rate. 


However, if the ‘Combustionflow’ model calculations were 





not limited by the hydrogen and oxygen flow rate 
limitations, then the maximum air flow rate becomes 80% of 


the choked air system flow rate. 
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6. Testing Time 
The maximum air flow rate computed above is sent to 


the ‘systemruntime’ model shown in figure (18). Here, the 





air supply system parameters are applied to the maximum air 





flow rat determined by the ‘Combustionflow’ model. 
Specifically, the air supply tank volume is evaluated to 
find the amount of time the maximum air flow rate can be 
supplied, while maintaining the tank pressure above the 
minimum value as specified in the user variables section of 


the program. The target final tank pressure is found by 





adding 250 psi to the stagnation pressure loop variable 
(unless this sum exceeds the supply air tank pressure). 


The model then returns the calculated run time to the main 





MATLAB program. 


At this point, the main program determines if the 


testing time is sufficient as discussed earlier in section 











T..B..2% If the run time does not meet the specification, 

















another SIMULINK model, ‘Makesystemruntime’ shown in figure 











(19), is called to lower the air mass flow rate until the 





minimum testing time requirements are met. 
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7. Facility Nozzle Throat Area and Area Ratio 





The final part of this program calculates the total 
mass flow rate and subsequently determines the facility 


nozzle parameters. The program sums the mass flow rates of 














the thr gases, Hz, O2 and Air to determine the total mass 


flow rate that the facility nozzle will experience. 


Since this is a choked condition at the facility 








nozzle throat, once again equation (9) is used. In this 
model, equation (9) is solved for throat area. The value 
is stored by the program in the ‘Astar’ matrix. The next 


step uses th ‘findmyexitmachnumber’ model shown in figure 





(20). This model uses a reference input for area ratio, 
and then applies Newton’s method to solve for the exit Mach 
number. When the matching Mach number is found the model 


stops and returns the exit Mach number to the main program. 








The process is then repeated for each of stagnation 
pressure and temperature conditions that were found to be 


valid in the atmosphere modeling process. 
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IV. RESULTS 


A. AIR PRESSURE SOURCE 
Several series of calculations were performed to 


determine the combustor inlet stagnation pressures. These 





calculations, depicted in figure (21), show that the 





maximum combustor stagnation pressure that occurs in the 
dynamic pressure flight window ranges roughly from 340 to 


2560 psi for the region of concern (flight Mach numbers 1 














to 10). 
x 10° Altitude vs Mach number 
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Figure 21. Region of Operation in the Flight Window 
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The graph above shows the black dynamic pressure lines 
forming a flight envelope dotted with red and blue 


asterisks. 


The upper black line represents a dynamic pressure of 





500 psf and the lower one is 2000 psf. The red “lines” are 

















lines of constant stagnation pressure. The blue lines ar 
lines of constant stagnation temperature. The upper red 
line is 340 psi and the lowest “broken” red line is about 





2560 psi combustor stagnation pressure. One aspect of the 








graph that should not be overlooked is that at a stagnation 
pressure (Pt) of 2563 psi, the full dynamic pressure 
(Paynamic) Window (from 500 to 2000 psf) can be replicated up 
to nearly Mach 7.5. Also the most restrictive condition of 
337.5 psi and 4400 °R and the mass flow rate is limited by a 
choked oxygen system flow. The corresponding air flow rate 
was determined to be about 17.1 lbm/s, which yields a test 


time of about 4 s. As discussed earlier, the minimum 





acceptable run time was set at 30 s. To achieve this 
minimum run time value without changing the piping systems, 


the air tanks would have to be 8 times larger than the 





current size if the air mass flow rate was kept at 17.1 





lbm/s. Conversely if the air supply volum wer kept 


constant, the air flow rate would have to be reduced to 1.9 





lbm/s to achieve 30 s test time. Reducing the mass flow 
rate by installing an orifice in the air supply line is 
certainly the cheapest solution in the short run but the 


reduction in system mass flow rate would require a smaller 











prototype combustor. A smaller combustor would introduce 





other problems that could lead to inaccurate test results. 





In order to minimize the cost and avoid some of the 


problems involved with very small scale modeling, a 
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compromise is proposed. If the volume of the air tanks is 


doubled, then the air system mass flow can be doubled to 





nearly 4 lbm/sec, and becaus th air mass flow rate at 





this temperature was only 65% by weight of the oxidizer, 





the total mass flow rate was about 6 lbm/sec. 


By increasing air storage volume, the RPCL would be 
able to perform hypersonic testing throughout the dynamic 


pressure flight window. 


B. VITIATOR DESIGN 


The design of a vitiated heater is without a doubt the 








most challenging part of this design. In the case of the 
RPCL design, the vitiator must be able to safely contain 
2560 psi and heat the air to 4400 °R. A cooled vitiator is 


necessary. One design is discussed by Hashimoto and 





Yoshida (Ref. 12). Their model was proven to operate at 
about 4 lbm/sec, 3000 °R and nearly 900 psi. While these 


conditions are not sufficient to support the conditions 





desired in this proposal, the temperature distribution 





profile was of interest. Specifically with a relatively 





small design they showed a 5-15% temperature stability ina 


short mixing distance. This was achieved by injecting the 





Oz and Hz through 8 injectors. They also designed the 
hydrogen and oxygen flow rates to be injected at different 


velocities, which results in shorter lengths to achieve 





macroscopic mixing. This success proves their design 


process and will aid in the design of a vitiator for the 





NPS hypersonic test stand. 
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Cc. FACILITY NOZZLE 

The facility nozzle design is one feature of the NPS 
test stand that will need to be adjustable, or several 
different nozzles will need to be available to support 
testing at different Mach numbers. This is due to the fact 
that the maximum total mass flow rate (not including test 
time restrictions) that RPCL can support was found to vary 
from about 3 to 45 lbm/sec, as shown in figure (22), over 


the entire dynamic flight window. 


Mass Flow Rate vs Stagnation Temperature and Stagnation Pressure 
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Figure 22. Mass Flow Rates versus Stagnation Temperatures at 
various Stagnation Pressures 


This resulted in correspondingly large variations in 
facility nozzle throat area, A’. Finally, the variation of 
A” at a constant mass flow rate was determined. The lowest 
system flow rate was used to normalize each ratio of mass 
flow rate to A’. Figure (23) shows that for the most 
restrictive condition, which corresponds to the lowest Pr 


and highest Tt, (337.5 psi and 4400 R), the facility nozzle 
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throat area varies by less than 0.03 in. This small 








variation will allow for the proposal of a facility nozzle 
design that consists of a common converging section with 
modular diverging extensions to allow for the supersonic 
area ratio to be adjusted to meet the exit Mach number 


requirements. 


Facility Nozzle Throat Area Variation 
at a constant mass flow rate at 2560 psi 
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Figure 23. Facility Nozzle Throat Area Variation 
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V. CONCLUSION 





In order to support a hypersonic blow-down facility at 


NPS’ RPCL, the air supply tank volume needs to be increased 


by at least 100% of its current capacity, (which is 160 ie 





at: 3000. psi; Ob Over’ 32,500: FE Of Ate ae. STP). This will 
allow testing of scramjet engines throughout a dynamic 
pressure flight envelope of 500 to 2000 psf in the Mach 
number range of 5 to 7.5. Additionally, the installed 


regulator will need to be bypassed to achieve stagnation 





pressures greater than 1000 psi. 


A method was developed and used to perform this 


evaluation. This method was implemented into a software 





suit developed to perform the calculations allowing the 





user to access multiple parameters and investigate the 


proposed hardware system changes to the design results. 








The preliminary requirements for the vitiator and 


facility nozzle were also found and basic solution concepts 





were proposed. A vitiated heater design burning Hz and Op, 


with evenly spaced multiple gas injection ports will 





generate the uniform temperature profile that is needed for 





this testing. Also since there are only small variations 








in the nozzle throat conditions, a modularized facility 





nozzle design is feasible, with diverging section 





extensions to allow the target exit Mach number to be 


adjusted. 


With minor support system modifications, the reality 


of a hypersonic test facility is within reach. 
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figure(1); 
title ('Altitude vs Mach number'); 
xlabel('Mach number'); 
ylabel('Altitude [ft]'); 

hold on 


for i=1:length(dynamp) ; 
dynampress=dynamp (i); 
j=1;Mach=[];alt=[]; 
for M=1:0.5:flightMmax; 
sim('atmosphericsolver'); 
if altout(end)<altmax; © Time limit of solver, if true then a 
solution was found 
Mach (j)=M; alt (3) =altout (end) ; 
j=jt+l; 





end 
end 
plot (Mach,alt,'k"'); 
end 


stagpress=stagplow; i=1; 
while stagpress<=Pairtank; 
for M=1:0.1:flightMmax; 

sim('Stagnationpressure'); 
dynampress=500; 
sim('atmosphericsolver'); 
highlimit=altout (end) ; 
dynampress=2000; 
sim('atmosphericsolver'); 
lowlimit=altout (end); 
if altoutp (end) <highlimitéaltoutp (end) >lowlimit; 





if i==1; 
stagp(i)=stagpress; 
i=it+l; 
end 
if i>lé&stagpress>stagp(i-1); % Records the stagnation 





pressures that fall within the defined flight window 
stagp(i)=stagpress; 
i=it+l; 
end 
plot (M,altoutp(end),'r*'); 
if machc (end) >machcmax; 


machcmax=machc (end) ; & For the flight window this 
specifies the maximum Combustor Mach number 
end 
end 
end 
stagpress=1.5*stagpress; 


end 


stagtemp=stagtlow;i=1; 
while stagtemp<=stagthigh; 
for Alti=0:2500:altmax; 
sim('Stagnationtemperaturebetter'); 
if mach(end)<flightMmax; % Solution found if true, limit of 
solver and valid combustor Mach number 
M=mach (end) ; 
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Records the stagnation 


pressures that fall within the defined flight window 


Q 
6 


For the flight window this 


, 


, 


° 
“6. 


a 


, 
, 


i a 
, 





1, 


stagtemp 
stagtemp 
machc (end) 


a 


, 


i 
R] 


[ 


2000 
sim('atmosphericsolver') 


lowlimit 


500 
sim('atmosphericsolver') 


highlimit 


=altout (end) 
altout (end) 


stagt (1) 
stagt (1) 
machcmax 
specifies the maximum Combustor Mach number 
7 


if i>l&stagtemp>stagt (i-1); 


plot (mach (end),Alti,'b*'); 
if machc (end) >machcmax 


if i==1 


end 
end 


if Alti<highlimitéAlti>lowlimit 
end 





dynampress 
dynampress 
end 

length (stagt 


1 





length (stagp) 


end 


stagp (x); 


=1 
for y 


end 
stagtemp=stagtemp+ (stagthigh-stagtlow) /10 


B>< >< ><> <>< ><> KOK OKO KOK OOK OKO KOKO KOKO KOOKS 
Facility subsystem parameters evaluation % 
B><>< >< >< ><> <> <> <> <DD> DKDOKDODOOKO KOKO KS 


end 
hold off 
for x 
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Diam(i); 


Dia 
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1, 


Press 
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Loop for the number of gas supply 
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fo) 


1, 


¥ 


gamma (i) 


gam 


1, 


0.8*mxmdot (end) 


length (GasConst) 


GasConst (i) 
sim('Chokedflowrate') 


mdotchoke (i) 
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for 2 
end 


systems 


oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 


oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 
oe 


63 


Calculates the O02 and H2 flow rates required to give the 


desired vitiator temperatures 
and returns the density of the mixture 











, 
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1, 


AR=1.0 


, 


Machfnout (end) 


Paral 
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the conditions in the flight window 
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Tt conditions. 
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